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Abstract 



We show the results of the fluctuation analysis applied 
to the off-source areas from the two 1 Million second Chan- 
dra Deep Fields, including our new results on the Chandra 
Deep Field-South (CDF-S) in the 0.5-2 keV band in addi- 
tion to those on the Chandra Deep Field-North (CDF-N), 
which have already been reported. The distribution of the 
X-ray counts in cells has been compared with the expec- 
tation from the Log N — Log S model to constrain the 
behavior of the source number density down to a factor 
of several lower than the source-detection limit. We show 
that our results are insensitive to the non-uniformity of 
the non X-ray background (NXB). Our results show that 
the number counts in the soft band (0.5-2 [keV]) continue 
to grow down to 5a, ~ 7 X 10 -18 [ergs -1 cm -2 ], possibly 
suggesting the emergence of a new population and they 
agree well with a prediction of star forming galaxies. 

Key words: galaxies: active — galaxies: evolution — (cos- 
mology:) diffuse radiation — X-rays: diffuse background 



1. Introduction 

The number densities of X-ray sources as a function of flux 
(the so called the Log N — Log S relation) is one of the 
key constraints for models of the X-ray source population. 
Most of the "Cosmic X-ray Background" (CXRB) inten- 



1993, Georgantopoulos ct al. 1993; Gendreau et al. 1998 



Yamashita et al. 2000; Pcrri fc Giommi. 2000). The source 



sity has now been resolved with "Chandra" ( Muhsotzky et 
al. 2000| ; |Rosati et al. 200l| ; |Brandt et al. 200l| ) to the ex 



tent that major uncertainties in the fraction of the CXRB 
which have been resolved into individual sources, i.e, how 
much of it remains to be explained, lie in the absolute in- 
tensity of the CXRB and field-to-field fluctuations due to 
cosmic variance. In terms of the origin of the CXRB, the 
faintest sources in the Chandra Deep Fields are becoming 
less interesting. However, constraints on number counts 
at the faintest possible fluxes provide a new view on the 
nature and evolution of the X-ray emitting sources. 

Fluctuation analysis is a strong tool for constraining 
the number densities of sources which are fainter than 
the source detection limit. This technique, which has been 



counts inferred from these analyses have turned out to 
be consistent with those from resolved sources in deeper 
observations later. Initial results of our work on CDF-N 



pioneered in radio astronomy (e.g. Condon 1974), has also 



been successfu lly applied to the X-ray data from previous 
missions (e.g. Hamilton fc Helfand 1987 ; Hasinger et al 



has been reported in Miyaji fc Griffiths (2002) (hereafter, 
MG02). In these proceedings, we also report our new result 
on the 0.5-2 keV result on CDF-S, which gave a consistent 
result, but more tightly constrained. 

We also discuss a number of technical issues which have 
not been explained in MG02. 

2. Image Simulation and Statistics 

The analysis has been made by comparing the count dis- 
tribution of 4" x 4" image cells (P(D) distribution) in 
the real data and the model. Because the point spread 
function (PSF) and exposure varies over the ACIS field of 
view (FOV), analytical modeling of the P(D) distribution 
cannot be used. Also because the probability distribution 
of the P(D) curve for a given model involves the fluctua- 
tions of both photon counts and source counts, it does not 
follow the simple Poisson statistics. Thanks to the rapid 
improvements in computing power in recent years, massive 
Monte-Carlo simulations can be made on our desktops for 
(1) modeling the P(D) distribution, and (2) searching for 
the 90% confidence range of the Log N — Log S relation. 

For this purpose, we have developed an image simu- 
lator (qimsim), which can quickly simulate the summed 
image from multiple pointings of Chandra observations 
and allows one to analyze the statistics of simulated im- 
ages. Qimsim generates randomly placed point sources for 
a given Log N — Log S model. For each source gener- 
ated at a sky position and for each pointing, it generates 
events based on the exposure map and the off-axis angle 
from the pointing direction. This is repeated for all the 
pointings used for the summed image. We assume that 
the background events are dominated by non X-ray (par- 
ticle) background (NXB). The NXB is generated using a 
separate "NXB distribution map" (see below) in such a 
way that the total of the simulated source and the NXB 
counts is equal to that of the real data. 

To assess the significance of the difference between the 
model P(D) (the mean of the P(D) curves from the 1000 
simulations) and that of the observation, we have used 
Monte-Carlo simulations. For this purpose, we can use any 
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C-Statistics 
(or any statistical measure of model-data deviation) 



high NXB and limited area of analysis. In these proceed- 
ings, we concentrate our discussion on the 0.5-2 keV band. 
Fig. U shows the 0.5-2 keV Log N — Log S relations for 
the CDF-N fluctuation, CDF-N resolved sources (MG02) 
and the new fluctuation result from CDF-S. [] The method 
of the CDF-S analysis is essentially the same as that of 
CDF-N reported in MG02. Since the pointing directions 
are much more concentrated in the CDF-S data than the 
CDF-N case, we used the single circle of 5 '2 around the 
exposure-weighted mean pointing direction for our analy- 
sis. We have also included the resolved source constraint 
at = 9 10 -17 [ergs~ 1 c m~ 2 ] from the resol ved source 



count and the la error from Rosati ct al. (2001) (see MG02 



for the detailed explanation on incorporating the resolved 
source constraint into analysis). Fig. ^ shows that the faint 
number counts of CDF-N and CDF-S are consistent with 
each other, but CDF-S gave a better constraint, probably 
because we can utilize a larger area for analysis. 



Figure 1. The assessment of acceptability of a Log N — 
Log S model using the Monte- Carlo simulations is illus- 
trated. The solid curve is the distribution o/C s i m (see text) 
for the 1000 simulations. The spread o/C s i m is caused by 
the statistical fluctuations of both the X-ray sources and 
the detected events. If C s i m exceeds C bs in 10% of the 
simulations or more (hatched area), we have accepted the 
model and included in drawing the 90% confidence fluctu- 
ation "horn'. 



quantity measuring the difference between the model and 
the data. For this purpose, we use the Caster modifica- 



tion of the Cash-C statistics ( Arnaud fc Dorman 2001 



Cash 1979 ; see also MG02). The acceptance probability 



has been assessed as as follows: 

1. Calculate the distribution of the 1000 C s i m values be- 
tween the model P(D) and the 1000 simulated P(D) 
distributions derived from the same set of model pa- 
rameters. 

2. Calculate C t, s between the data and the model. 

3. The fraction of the C s ; m values which exceed C b s is 
taken as the model acceptance probability. 

4. If this fraction is > 0.1, we accept it. From the points in 
the model parameter space which have been accepted, 
we calculate the the 90% confidence "horn" . 

This procedure is illustrated in Fig. [l]. 

We plan to use the simulator for further studies involv- 
ing angular correlation functions. 

3. Fluctuation Results, AGNs and Galaxies 
3.1. Results 

As we see in MG02, we could only obtain loose constraints 
on the 2-10 keV band for CDF-N, mainly because of the 
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Figure 2. The derived Log N — Log S relations ( 90% con- 
fidence "horns") in the soft band in the CDF North and 
South are compared with AGN and galaxy number count 
models. 

Predictions from the population systhesis model (B) by 
Gilli et al. (20Q~T\ for the total AGN population and the 
obscured (N-r > 10 22 [cm~ 2 ]). Predicted number counts of 
galaxies using the cosmic star-formation history and two 
models ( Gaussian and Peak-M) of evolution of X-ray bi- 
naries by Ptak et al. (200~T\ are also plotted. If mode B 



of Gilli et al. (2001 i represents the correct behavior of ab- 
sorbed and unabsorbed AGNs, the fluctuation constraints 
suggest the emergence of an additional population, proba- 
bly from these galaxies. 



1 Machine readable ASCII tables of the fluctuation 
ho rns for our analysis on CDF-N /S can be found at 



http://astrophysics.phys.cmu.edu/~deepxray/tables 
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We have overplotted the 0.5-2 keV Log N — Log S pre- 
diction from the AGN population synthesis model com- 
posed of unabsorbed and absorbed AGNs based on by 
Gilli et al. (2001)| (model B). The plotted model have 



been slightly modified from the original as described in 
Rosati et al. (200f). The curves for the total and ab- 
sorbed (Ah > 22 [cm -2 ]) AGN populations are plotted 
Fig .|2| shows that the AGN counts drop below S x ~ 10 
[ergs -1 cm -2 ] even if the emergence of the obscured AGN 
population is taken into account, while the source counts 
from the fluctuation analysis continue to grow. If their 
model B closely represents the true behavior of the SXLF 
of the AGN population, we are probably seeing the emer- 
gence of a new population of faint X-ray sources. This ex- 
cess may be contributed by AGNs with low intrinsic lumi- 
nosities (at L x 10 41 ' 5 [ergs -1 ]) and/or X-rays from star- 
formation activities through supernova remnants and low- 
/high- mass X-ray binaries. In view of this, we have also 
overplotted a model prediction of X-ray number counts 



based on the cosmic star-formation rate Ptak et al. (2001) 
for the two models of binary evolution (Gaussian and 
Peak-M) from phosh fc White (200 1)| . Within the uncer- 
tainties in the AGN source counts and model construction, 
these predictions gave approximately the number counts 
from our fluctuation analysis. 

3.2. Sensitivity to Non-Uniform NXB 
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Figure 3. Sensitivity to the NXB non uniformity has been 
assessed. The solid curve assumes that the NXB is uni- 
form over the active region of the CCD and dotted curve 
assumes that the NXB has a 20% gradient across the 
ACIS-I CCD chips. The comparison of these two shows 
that our analysis is not sensitive to the non-uniform NXB 
level. 



It has been reported that the non X-ray background 
(NXB) varies over the ACIS CCD-chip at the level of 
- 20 - 30% across the ACIS-I chip for E > 5 keV, while 
the variation is much less for lower energies. []. While we 
treated NXB uniform over the active region in our anal- 
ysis, one may worry that this NXB non uniformity could 
cause an overestimate of the derived Log N — Log S rela- 
tion. Thus we have also made the analysis assuming that 
the NXB has a 20% gradient along the CHIP-Y coordi- 
nate by adjusting the NXB distribution map. The results 
for the uniform and "slope" NXB are compared in Fig. 
||. Fig. H shows that the systematic error caused by this 
effect is much smaller than statistical errors. 



4. Summary 

Fluctuation analyses have been made on the two 1 Ms 
Chandra observations on the CDF-North and South in 
order to constrain the behavior of the Log N — Log S 
at fluxes fainter than the source detection limit. We 
also show our new results on CDF-S for the soft band. 
We have used Monte-Carlo simulations to obtain model- 
predicted P{D) curves and evaluate the 90% confi- 
dence constraints. 

The number counts in the soft band (0.5-2 [keV]) con- 
tinues to grow down to S x ~ 7 x 10 -18 [ergs -1 cm -2 ], 
suggesting the emergence of a new population, possibly 
X-rays from star-formation activity. 
The impact of non uniformity of the non X-ray back- 
ground on our Log N — Log S is much smaller than 
the statistical errors. 
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